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Abstract 

The OPERA neutrino detector in the underground Gran Sasso Laboratory 
(LNGS) has been designed to perform the first detection of neutrino oscilla- 
tions in direct appearance mode through the study of the — > i/^ channel. 
The hybrid apparatus consists of an emulsion/lead target complemented by 
electronic detectors and it is placed in the high energy long-baseline CERN 
to LNGS beam (CNGS) 730 km away from the neutrino source. Runs with 
CNGS neutrinos were successfully carried out in 2008 and 2009. After a 
brief description of the beam, the experimental setup and the procedures 
used for the analysis of the neutrino events, we describe the topology and 
kinematics of a first candidate Uj. charged-current event satisfying the kine- 
matical selection criteria. The background calculations and their cross-check 
are explained in detail and the significance of the event is assessed. 

Keywords: 
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1. Introduction 



Neutrino oscillations were anticipated nearly 50 years ago[H but were 
observed much later [2j. Several experiments carried out in the last decades 
with atmospheric and accelerator neutrinos, as well as with solar and reactor 
neutrinos, have contributed to our present understanding of neutrino mixing 
(see e.g. j3| for a review). 

Atmospheric neutrino oscillations, in particular, have been studied by the 
Kamiokande 8, Super-Kamiokande Q, MACRO and S0UDAN2 Q ex- 
periments. The CHOOZ [3] and Palo Verde jsf reactor experiments excluded 
indirectly the i/^ channel as the dominant process in the atmospheric 

sector. Two accelerator experiments have already confirmed the oscillation 
hypothesis in disappearance mode, K2K [9|] and MINOS 10|. However, the 
direct observation of flavour transition through the detection of the corre- 
sponding lepton has never been observed. Appearance of will prove un- 
ambiguously that — ?■ oscillation is the dominant transition channel at 
the atmospheric scale. This is the main goal of the OPERA experiment 



11, 12, 13 



14 neutrino beam 



with the high energy long baseline CNGS 
from CERN to the Gran Sasso Underground Laboratory of INFN. The ob- 
servation of interactions of CNGS neutrinos in the OPERA detector have 



already been presented [151 . 

In this paper, we report the observation of a first r-lepton candidate 
event, possibly produced by an oscillated z/^, consistent with the kinematical 



selection criteria [12|, ll3| . We evaluate the expected numbers of signal and 
background events in the data sample analysed so far, the probability that 
the event is due to a background fiuctuation and the statistical significance 
of the observation. 



2. The OPERA detector and the CNGS beam 

The challenge of the OPERA experiment is the detection of the short- 
lived T lepton (cr = 87 /im) produced in the charged-current (CC) interaction 
of Vr- This sets two requirements difficult to conciliate: a large target mass 
to collect enough statistics and a very high spatial accuracy to observe the r 
lepton, the decay length of which being in the millimetre range at the CNGS 
beam energy. In OPERA, neutrinos interact in a large mass target made of 
lead plates interspaced with nuclear emulsion films acting as high-accuracy 
tracking devices. This kind of detector is historically called Emulsion Cloud 
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Chamber (ECC). It was successfull y u sed to establish the first evidence for 



charm in cosmic rays interactions |l6l| and in the DONUT experiment for 



the first direct observation of the Ur [13]. OPERA is an hybrid detector [18| 



made of a veto plane followed by two identical Super Modules (SM) each 
consisting of a target section of about 625 tons made of 75000 emulsion/lead 
ECC modules or "bricks", of a scintillator Target Tracker detector (TT) 
to trigger the read-out and localize neutrino interactions within the target, 
and of a muon spectrometer. A target brick consists of 56 lead plates of 1 
mm thickness interleaved with 57 emulsion films and weighs 8.3 kg. Their 
thickness along the beam direction corresponds to about 10 radiation lengths. 
In order to reduce the emulsion scanning load. Changeable Sheets (CS) film 
interfaces have been used. They consist in tightly packed doublets of emulsion 
films glued to the downstream face of each brick. 

Charged particles from a neutrino interaction in a brick cross the CS and 
produce signals in the TT that allow the corresponding brick to be identified 
and extracted by an automated system. Large ancillary facilities are used 
to bring bricks from the target up to the automatic scanning microscopes 



at LNGS and various laboratories in Europe and Japan [19|, |20|. Extensive 



information on the OPERA detector and ancillary facilities are given in 18|, 



21 



OPERA is exposed to the long-baseline CNGS z-'^ beam [ij] from CERN 
to LNGS, 730 km away. The beam is optimized for the observation of Vr CC 
interactions. The average neutrino energy is ~17 GeV. The contamination 
is 2.1% in terms of interactions; the Pe and contaminations are lower than 
1%, while the number of prompt Vr is negligible. With a total CNGS beam 
intensity of 22.5 10^^ protons on target (p.o.t.), about 2s300 neutrino events 
would be collected. The experiment should observe about 10 Vr CC events 
for the present ^m\^ allowed region with a background of less than one event. 
In 2006 and 2007, first interactions were recorded during the commissioning 
phase of the CNGS. Significant amount of data were collected starting in 2008 
[15]. We report here on results based on part of the data taken during 2008 
and 2009 runs, which amounted respectively to 1.78 and 3.52 x 10^^ p.o.t., 
respectively. Among the 10122 and 21428 on-time events taken in 2008 and 
2009, respectively 1698 and 3693 events were identified as candidate neutrino 
interactions in the target bricks. 
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3. Event selection and analysis 



The appearance of the r lepton is identified in OPERA by the detection 
of its characteristic decay topologies, either in one prong (electron, muon or 
hadron) or in three pron gs. 

A software algorithm 22] selects CC and neutral-current (NC) events in- 
side the target. When a trigger in the electronic detectors is compatible with 
an interaction inside a brick, a software reconstruction program processes the 
electronic detector data to select the brick with the highest probability to 
contain the neutrino interaction vertex. The contamination by interactions 
occurring outside the target is evaluated by Monte Carlo simulations to be 
4.8%. Bricks designated by the finding algorithm are removed from the tar- 
get wall and exposed to X-rays for CS-to-brick alignment and the CS doublet 
is detached from the brick and developed underground. The corresponding 
CS films are then analysed to validate the designated brick. The information 
of the CS is later used for a precise prediction of the position of the tracks 
in the most downstream films of the brick, helping the tracks scan-back in 
the subsequent vertex finding procedure. Validation of the selected brick 
is achieved if at least one track compatible with hits reconstructed in the 
electronic detectors is detected in the CS films. The efficiency of the brick 
finding procedure from the analysis of the 2008 run is found to be 77% after 
the most probable or, when required, the two most probable bricks have been 
analysed. This efficiency reaches 83% in a subsample where up to 4 bricks 
per event were processed. 

After a brick has been validated, it is brought to the surface to be ex- 
posed to high-energy cosmic-rays for a precise film-to-film alignment. The 
brick emulsion films are then developed and dispatched to the various scan- 
ning laboratories. All tracks measured in the CS are sought in the most 
downstream films of the brick and followed back until they are not found in 
three consecutive films. The stopping point is considered as the signature 
either for a primary or a secondary vertex. The vertex is then confirmed 
by scanning a volume with a transverse size of 1 cm^ in 15 films in total, 5 
upstream and 10 downstream of the stopping point. The present overall loca- 
tion efficiency averaged over NC and CC events, from the electronic detector 
predictions down to the vertex confirmation, amounts to about 60%. 

A further phase of analysis is applied to located vertices, called decay 
search procedure, to detect possible decay or interaction topologies on tracks 
attached to the primary vertex. It is also searched for extra tracks from 
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neutral decays, interactions and 7-ray conversions. When secondary vertices 
are found in the event, a kinematical analysis is performed, using particle 
angles and momenta measured in the emulsion films. 

For charged particles up to about 6 GeV/c, momenta can be deter- 
mined using the angular deviations produced by Multiple Coulomb Scatter- 
ing (MCS) of tracks in the lead plates [23.!]. This method gives a momentum 
resolution better than 22% for charged particles with momenta lower than 6 
GeV/c, passing through an entire brick at normal incidence. For higher mo- 
mentum particles, the measurement is based on the position deviations. The 
resolution is better than 33% on 1/p up to 12 GeV/c for particles passing 
through an entire brick at normal incidence. Momenta of muons reaching 
the spectrometer are measured with a resolution of 20% up to 30 GeV/c and 
the sign of their charge is determined jj^ . 

The 7-ray search is performed in the whole scanned volume by checking 
all tracks having an impact parameter (IP) with respect to the primary or 
secondary vertices lower than 800 fim. The angular acceptance is ±500 mrad. 
The 7-ray energy is estimated by a Neural Network algorithm that uses the 
combination of the number of segments and the shape of the electromagnetic 
shower and also the multiple Coulomb scattering of the leading tracks. The 
same method is applicable to electrons from CC interactions. The 7- 
ray attachment to a vertex is based on the angular resolution of the first 
measured segment of the e~^e~ pair. The precision in the 7-ray attachment 
to a vertex depends on the conversion point inside the 1 mm thick lead plate. 
The typical resolution for 7-ray angle determination is 5 (3) mrad at 1 (3) 
GeV. 

The detection of decay topologies is triggered by the observation of a 
track with a large impact parameter with respect to the primary vertex. The 
IP of primary tracks is smaller than 12 fim after excluding tracks produced 
by very soft particles. 

The results presented in this paper come from the decay search analysis of 
a sample of 1088 events of which 901 are classified as CC interactions. This 
sample, taking into account the event location efficiencies, corresponds to 
1813 interactions occurring in the target, namely about 35% of the 2008 and 
2009 data sample, or, rescaled to the beam integrated intensity, to 1.89x10^^ 
p.o.t.. 

Charmed particles have similar lifetimes as r leptons and, if charged, 
share the same decay topologies. Being produced in i/^ CC interactions, 
the location of their primary vertex is easier than that of muon-less Uj. CC 
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interactions where the r lepton decays into hadrons or into an electron. Still, 
the primary vertex being located, the finding efficiency of the decay vertices 
is similar for both types of particles, although the selection criteria are more 
stringent in the case of the r analysis that is aiming at minimal background. 
In the sample of z/^ CC interactions that have been searched for, a total 
of 20 charm decays have been observed that survived selection cut-offs, in 
agreement with the number expected from a MC simulation, 16.0±2.9. This 
demonstrates that the efficiency of the search for short-lived decay topologies 
is understood. Among them, 3 have a one-prong topology where 0.8 ± 0.2 
are expected. The background in the total charm events sample is about 2 
events. 



4. Candidate event topology and track kinematics 

The decay search procedure applied to the 2008-2009 data sample anal- 
ysed so far yielded one candidate event with measured characteristics fulfill- 
ing the selection criteria a priori defined for the v-r interaction search. The 



cuts used in this selection are those defined in the experiment proposal [12 
and its addendum 1 



For decays of the r to a single charged hadron, it is required at the 
secondary vertex, when not in the primary vertex lead plate, that: 

• the kink angle ^fcmfc is larger than 20 mrad; 

• the secondary vertex is within the two lead plates downstream of the 
primary vertex; 

• the momentum of the charged secondary particles is larger than 2 
GeV/c; 

• the total transverse momentum {Pt) of the decay products is larger 
than 0.6 GeV/c if there are no photons emitted at the decay vertex, 
and 0.3 GeV/c otherwise. 

At the primary vertex, the selection criteria are: 

• there are no tracks compatible with that of a muon or an electron; 

• the missing transverse momentum (P™**'^) is smaller than 1 GeV/c; 
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• the angle $ in the transverse plane between the r candidate track and 
the hadronic shower direction is larger than 7r/2. 

The muon-less neutrino interaction reported in detail here occurred in a 
brick situated in wall 11 of the first SM and well inside the target: 3 bricks 
from the top and 24 bricks from the left. This allows a comprehensive study of 
the event: tracks can be followed over large distances and secondary vertices 
including electromagnetic showers searched for in a large volume. The event 
has been independently measured with the European and Japanese scanning 
systems. Measurements are consistent and their averages are considered in 
the following. 

The primary neutrino interaction consists of 7 tracks of which one exhibits 
a visible kink. Two electromagnetic showers caused by 7-rays have been 
located that are associated with the event. Fig. [1] shows a display of the 
event. 

• Track 1 exits from the primary interaction brick and is not found in the 
brick immediately downstream. It is left by a particle of momentum 
(0.78lQ;j^Q) GeV/c most likely interacting in the target tracker between 
both walls; 

• Track 2 is left by a heavily ionizing particle. From its residual range 
(32.0 ±0.5) g cm-2 and the value of j9/3 = (0.321°;?}) GeV/c measured 
on the upstream half of the track, the particle is identified as a proton, 
the kaon hypothesis being rejected with a CL. of 97%. The proton 
momentum resulting from the residual range corresponds to (0.60 ± 
0.05) GeV/c; 

• Track 3 is left by a particle which generates a two-prong interaction 4 
bricks downstream of the primary vertex. Its momentum is equal to 
(1.97l°:i) GeV/c; 

• Track 5 has been followed in wall 12 and disappears in wall 13 after a 
total distance shorter than 174 g cm~^. The particle has a momentum 
of il.30t'oil) GeV/c; 

• Track 6 stops in the primary brick. The range, ionisation and multiple 
scattering are that of a pion of very low momentum (0.36tQ;o9) GeV/c; 
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Figure 1: Display of the r candidate event. Top left: view transverse to the neutrino 
direction. Top right: same view zoomed on the vertices. Bottom: longitudinal view. 
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• Track 7 is not directly attached to the primary vertex and points to 
it with an IP = (43I43) A*^- starting point is separated from this 
vertex by 2 lead plates. Its origin is likely to be a prompt neutral 
particle. In the analysis, the momentum (0.49lg^3) GeV/c has been 
added to the total momentum at the primary vertex; 

• Track 4 exhibits a kink topology with an angle of (41±2) mrad after a 
path length of (1335±35) /im corresponding to 15 cr would it be that of 
a r lepton. The expected 7 factor from the kink angle is approximately 
25. It is the parent track of a secondary interaction or decay. Both the 
kink angle and the path length satisfy the selection criteria (Figs 
andEJo); 

• Track 8, the kink daughter track, is left by a particle of a high mo- 
mentum of (I2I3) GeV/c well above the 2 GeV/c selection cut-off and 
which generates a 2-prong interaction 7 walls downstream its emission 
vertex (Fig. Its IP with respect to the primary vertex is (55 ± 4) 
fim. All the tracks directly attached to the primary vertex match the 
vertex point within 7 //m. 

None of the tracks of the charged particles emitted at either vertex is 
compatible with being that of an electron. Based on momentum-range con- 
sistency checks, the probability that tracks 1, 5 and 6 are left by muons is 
estimated to less than 0.001. The remaining tracks are identified as left by 
hadrons through their interaction. The probability of a soft muon with a 
slope larger than 1, which is undetectable in the brick and also in the TT, is 
1.5%. 

7-ray 1 has an energy of (5.6 ± 1.0(sif:at.) ±1.7{syst.)) GeV. The distance 
between its conversion point and the secondary vertex is 2.2 mm and the 
shower is compatible with pointing to it with a probability of 32%, the impact 
parameter being (7.5 ±4.3) ^m. Its probability to be attached to the primary 
vertex is less than 10~^, the IP being (45.0 ± 7.7) /im. 

7-ray 2 has an energy of {1.2±0. A{stat.)±OA{syst.)) GeV. It is compatible 
with pointing to either vertex, with a significantly larger probability of 82% 
at the secondary vertex, the IP being (22I22) f^^^ compared to 10% and 
(85±38) /im at the primary vertex. Its distance to both vertices is about 13 
mm. 
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Figure 2: MC distribution of: a) the kink angle for r decays, b) the path length of the 
T. c) the momentum of the decay daughter, d) the total transverse momentum Pt of the 
detected daughter particles of r decays with respect to the parent track. The red band 
shows the 68.3% domain of values allowed for the candidate event and the dark red line the 
most probable value. The dark shaded area represents the excluded region corresponding 
to the a priori tau selection cuts. 
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5. Kinematical analysis of the candidate event 

In the most probable hypothesis 7-ray 2 is emitted at the secondary ver- 
tex. The total transverse momentum Pt of the daughter particles with re- 
spect to the parent track is (0.47tai2) GeV/c, above the lower selection 
cut-off at 0.3 GeV/c. In the low probability hypothesis where 7-ray 2 is at- 
tached to the primary vertex, the systematic shift in Pt does not exceed 50 
MeV/c. This small effect is included in Fig. |2]i. 

The missing transverse momentum P™** at the primary vertex is (O.drtofr) 
GeV/c. This is lower than the upper selection cut-off at 1 GeV/c (Fig. [3^). 
The angle $ between the parent track and the rest of the hadronic shower 
in the transverse plane is equal to (3.01 ± 0.03) rad, largely above the lower 
selection cut-off fixed at 7r/2 (Fig. [3]d). The total hadronic momentum at the 
primary vertex, not including that of the parent assumed to be a r lepton, is 
about 5.5 GeV/c. The sum of the modulus of the momenta at the primary 
vertex is (24.3^32) GeV/c. This compares well with the total energy mea- 
sured in the TT by calorimetry, (44±12) GeV, owing for the fraction of the 
momenta taken by the neutrals invisible in the bricks. 

The invariant mass of 7-rays 1 and 2 is (120 ± 20{stat.) ± 35{syst.)) 
MeV/c^, supporting the hypothesis that they originate from a 7r° decay. 
Similarly the invariant mass of the charged decay product assumed to be 
a 7r~ and of the two 7-rays amounts to {640^11^ {stat.)'^lQ^{syst.)) MeV/c^, 
which is compatible with the p{770) mass. The branching ratio of the decay 
mode T -> p^i^r is about 25%. 

6. Background estimation and statistical significance 

The secondary vertex is compatible with the decay of a t~ into /i^(n7r°)z/^. 
The two main sources of background to this channel where a similar final state 
may be produced are: 

• the decays to a single charged hadron of charged charmed particles pro- 
duced in z/^ CC interactions where the primary muon is not identified 
as well as the c — c pair production in z/^ NC interactions where one 
charm particle is not identified and the other decays to a one-prong 
hadron channel; 

• the one-prong inelastic interactions of primary hadrons produced in i/^ 
CC interactions where the primary muon is not identified or in NC 
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Missing (GeV/c) (raj) 

Figure 3: Monte-Carlo distribution of a) the missing Pt at the primary vertex and b) the 
angle $ between the parent track and the rest of the hadronic shower in the transverse 
plane at the primary vertex for Vr CC interactions in black and for v^^ NC background in 
blue. The red band shows the 68.3% domain of values allowed for the candidate event and 
the dark red line the most probable value. The dark shaded area represents the excluded 
region corresponding to the a priori tau selection cuts. 



interactions and in which no nuclear fragment can be associated with 
the secondary interaction. 

The charm background produced in interactions in the analysed sam- 



ple amounts to 0.007±0.004 (syst.) events [l2|, ll3|]. That produced in 
interactions is less than 10~^ events. Current estimates of charm background 
are conservative since they are based on the original scanning strategy of the 
proposal. They do not include the additional reduction obtained by following 
all tracks up to their end points as was done for this candidate event. An 
improved evaluation of the charm background that exploits such additional 
information is in progress. 

A search was performed for hadronic activity seeking for highly ionising 
tracks associated with nuclear fragments and pointing to the secondary ver- 
tex. No such track is observed while the number of a-ray tracks in the vicinity 
of the event is found to be compatible with that expected from the lead natu- 



ral radioactivity 21[. The dominant background from hadron re- interactions 



has been evaluated with a FLUKA based Monte Carlo code [25[. About 160 
millions of events were simulated where 0.5-15 GeV tt"'", tt", , K~ and p 
are impinging on 1 mm of lead. The probability for a background interaction 
to occur over 2 mm of lead, the maximal decay length considered, and to 
satisfy the selection criteria for the reconstruction of the kink decay topology 
and its kinematics is (1.9 ± 0.1) x 10""^ per NC event. This probability de- 
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creases to (3.8 ±0.2) x 10~^ per NC event when taking into account the cuts 
on the event global kinematics. This leads to a total of 0.011 ± 0.006{syst.) 
background events when misclassified CC events are included. We assume 
a 50% error on the result of the FLUKA simulation in order to take into 
account possible systematic effects in the hadronic interactions modeling. 

The total background in the decay channel to a single charged hadron is 
therefore 0.018 ± 0. 007 {syst.). This new estimation based on an upgraded 
and more reliable simulation tool is 2.2 times larger than the value of the 



experiment proposal [12|, [13 



Three experimental cross-checks of the one-prong hadron background cal- 
culation have been performed, though so-far with small statistics. The tracks 
of hadrons from a sub-sample of NC and CC interactions have been followed 
far from the primary vertex over a total length of 8.6 m and no background- 
like interaction has been found in the signal region. This corresponds to a 
probability of occurrence over 2 mm of lead smaller than 1.54 x 10~^ at 90% 
C.L. Outside the signal region, a total of 7 single-prong interactions have 
been located with a kink angle > 20 mrad and a daughter transverse mo- 
mentum Pt > 0.2 GeV/c while 8 are expected from the simulation. Within 
the low statistics, both Pt distributions agree well. Finally, a brick has been 
exposed to a 4 GeV/c tt- beam and 119 interactions have been analysed so 
far. The forward and backward track multiplicities, the kink angle, the sec- 
ondary particles momentum and their Pj- have been checked against FLUKA. 
They agree within statistical uncertainties both in normalisation to ~10% 
and shape. (7.4 ± 1.1)% are one-prong interactions with a kink angle larger 
than 20 mrad where (7.6 ±0.3)% are predicted by the simulation. 

Finally, the expected number of events due to prompt ly^ background from 
Ds decays is 10~^ (10~^ per CC event). 

The probability that the expected 0.018 ±0.007(syst.) background events 
to the h~{mr^)h'r decay channel of the r fluctuate to one event is 1.8% 
(2.36 a). As the search for r- decays is extended to all four channels, 
f^'T^ri^fi, e'Urh'e, {mi^)vr and tt^tt^tt+z/t-, the total background then be- 
comes 0.045 ± 0.023(syst.). The probability that this expected background 
to all searched decay channels of the r fluctuate to one event is 4.5% (2.01 
a). 

At Am23 = 2.5 X 10 ^ eV^ and full mixing, the expected number of 
observed r events with the present analysed statistics is 0.54 ± 0.1?>{syst.) 
of which 0.16 ± 0.0'i{syst.) in the one-prong hadron topology, in agreement 
with the observation of one event. 
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7. Conclusions 



We reported the observation of a first candidate f-r CC interaction in the 
OPERA detector at LNGS. It was identified in a sample of events corre- 
sponding to 1.89 X 10^^ p.o.t. in the CERN CNGS beam. The assumed 
T~ lepton decays into h~(mT^)i'r and passes the selection criteria [l2[ 13|. 

The observation of one possible tau candidate in the decay channel h~{n^)h'r 
has a significance of 2.36 a of not being a background fluctuation. If one 
considers all decay modes included in the search the significance of the ob- 
servation becomes 2.01 a. The expected number of u^- events detected in the 
analysed sample is 0.54 ± 0.13{syst.). 

We therefore claim the observation of a first candidate CC interaction. 
Its significance, based on our best conservative knowledge of the background, 
exceeds two a. This does not allow yet claiming the observation of z/^ — )■ Ur 
oscillation. Given its sensitivity, the OPERA experiment will require the 
detection of a few more candidate events in order to firmly establish neutrino 
oscillations in direct appearance mode through the identification of the final 
charged lepton. 
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